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Abstract 

The electromagnetic properties of baryons containing a heavy quark are calculated at next-to- 
leading order in partially quenched heavy hadron chiral perturbation theory. Calculations are 
performed for three light flavors in the isospin limit and additionally for two light non-degenerate 
flavors. We use partially-quenched charge matrices that are easy to implement on the lattice. 
The results presented are necessary for the light quark mass extrapolation and zero-momentum 
extrapolation of lattice QCD and partially quenched lattice QCD calculations of heavy hadron 
electromagnetic properties. Additionally relations between the sextet electromagnetic form factors 
and transition form factors are derived. 
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I. INTRODUCTION 



Calculating physical properties of hadrons from QCD is a formidable challenge. The 
strong interaction dynamics involved is highly non-perturbative and results in the confine- 
ment of quarks and gluons into color-neutral hadronic states. Electromagnetic probes have 
often been sought experimentally as a means to glimpse the charge and magnetism distribu- 
tions of quarks within hadrons. Electromagnetic observables, such as the charge radii and 
electromagnetic moments, provide a physically intuitive glimpse at the structure of hadrons. 
For baryons containing 6 or c quarks, however, the experimental difficulties in measuring 
these observables are substantial. 

The dynamics underlying singly heavy baryons involves a rich interplay of opposite limits 
of QCD: heavy-quark physics and light-quark physics. Lattice QCD provides a first prin- 
ciples method for the calculation of QCD observables. These numerical simulations can 
determine the electromagnetic properties of baryons containing a heavy quark and thereby 
provide a physical picture of these baryons in terms of their quark degrees of freedom. Calcu- 
lation of singly heavy baryon observables (masses and weak transitions) have been pursued 
in quenched QCD (QQCD) QllSliSBQIi- Until now 0, partially quenched QCD 
(PQQCD) calculations have been very limited for any baryons. The appearance of quenched 
and partially quenched approximations to lattice QCD arises from numerically calculating 
the fermionic determinant. In QQCD, the time-costly numerical computation is radically 
simplified by replacing the determinant by a constant, whereas in PQQCD the determinant 
is calculated, but using larger quark masses for the quarks not coupled to external sources. 
In PQQCD, sea quark contributions are thus retained; while in QQCD, the sea quarks have 
been effectively discarded. Due to restrictive computational time, the valence quark masses, 
too, cannot be set to their physical values and one must extrapolate down from the values 
used on the lattice. 

To perform quark mass extrapolations, one must understand how QCD observables be- 
have as the quark masses vary. This can be done using chiral perturbation theory (xPT), 
which is a model- independent effective theory for low-energy QCD. This effective theory is 
written in terms of the pseudo-Goldstone bosons appearing from chiral symmetry break- 
ing ^ 11 1 . While these bosons are not massless, they remain light compared to the scale of 



chiral symmetry breaking and dominate low-energy hadronic interactions. The light quark 
masses appear as parameters of xPT and enable the determination of the quark mass de- 
pendence of QCD observables by matching onto the effective theory. Without external input 
(e.g. experimental data or lattice QCD data), the effective theory has little predictive power 
because symmetries constrain the types of operators in the Lagrangian, but not the val- 
ues of their coefficients. While providing a guide to extrapolate lattice QCD data, xPT 
can also benefit in turn from the determination of its low-energy constants (LECs) from 
these data. Heavy quark symmetry furthermore can be combined with chir al sy mmetry to 
describe hadrons formed from light quarks and a single heavy quark 

El HQ- 

For QQCD lattice simulations, quenched chiral perturbation theory (QvPT) has been 
developed to aid in the extrapolation in valence quark masses [isl EgI. llTlllilligt. There is, 
however, no general connection of QQCD observables to QCD because, for example, QQCD 
lacks an axial anomaly. Consequently QxPT contains operators that do not have analogues 
in xPT. Furthermore, the LECs of operators that have analogues are numerically unrelated 
to their counterparts in xPT. The sickness of the quenched approximation can be treated 
by utilizing partially quenched lattice simulations. For such simulations, partially quenched 
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chiral perturbation theory PQ^PT has been constructed to perform the extrapolation in 
both sea and valence quark masses 0, HHIil H Q . Unlike QQCD, PQQCD retains 
an axial anomaly and the flavor singlet field can be integrated out. In further distinction to 
QxPT, the LECs of xPT appear in PQ^PT, and hence PQQCD lattice simulations can be 
used to determine yPT parameters. Much work has resulted from this possibility, e.g 



m 



the baryonic sector see 0, H 111 111, H H H H H . 

In this work we determine the electromagnetic properties of baryons containing a heavy 
quark at next-to leading order (NLO) in partially quenched heavy hadron xPT. We de- 
termine the electromagnetic properties in three-flavor and two-flavor partially quenched 
simulations and use the isospin limit in the former. Such expressions allow one to perform 
the required extrapolation in light-quark masses. Additionally the electromagnetic LECs 
that appear in heavy hadron xPT can be extracted at NLO from lattice data and we ex- 
plore how the arbitrariness of the light-quark charge matrix can be exploited to determine 
the LECs. Furthermore we derive relations between the electromagnetic form factors and 
electromagnetic transition form factors of the sextet baryons. These relations follow from 
the decoupling of the heavy and light spin degrees of freedom. 

The organization of the paper is as follows. First in Sec. m we review the inclusion 
of baryons containing a heavy quark into PQxPT. Next, the charge radii and magnetic 
moments of the se = baryons are calculated in Sec. IIIIl The charge radii, magnetic 
moments and electric quadrupole moments of the S£ = 1 baryons are calculated in Sec. IIVI 
The corresponding results for SU{4\2) with non-degenerate light quarks are presented in 
Appendix^ For reference, in Appendix^ we give the electromagnetic properties of these 
baryons in SU{3) and SU{2) xPT. In Appendix O we use heavy quark symmetry to derive 
relations between the sextet baryon form factors; and in Appendix El we present results 
for and comment on the commonly used form of the partially-quenched charge matrix. The 
summary. Sec. |3 highlights the goal of understanding the electromagnetic properties of 
the triplet and sextet baryons from lattice QCD, and contrasts the situation with that of 
baryons formed from three light quarks. 



II. PQxPT FOR HEAVY HADRONS 

In this section, we present the formulation of PQxPT for heavy hadrons in SU{6\3). We 
begin by reviewing the set up of PQQCD. Next we describe the pseudo-Goldstone mesons 
of PQxPT, and finally include the heavy hadrons into this partially quenched theory. 



A. PQQCD 

The light-quark sector Lagrangian in PQQCD is given by 

9 

£= J2qj{z]p-m,).,q,. (1) 

j,k=i 

This differs from the QCD Lagrangian by the addition of six extra quarks; three bosonic 
ghost quarks, {u,d,s), and three fermionic sea quarks, {j,l,r), in addition to the light 
physical quarks {u, d, s). The nine quark fields transform in the fundamental representation 
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of the graded group SU(6\3) j35|, |36[. They appear in the nine- component vector 

q = {u,d,s,j,l,r,u,d,s)^. (2) 
These quark fields obey graded equal-time commutation relations 

?r(x)gf (y) - (-i)'''''^gf (y)gr(x) = ^'^^..^^(x - y), (3) 

where a, (3 and i,j are spin and fiavor indices, respectively. The graded equal-time commu- 
tation relations which vanish can be written similarly. The different statistics of the PQQCD 
quark fields are refiecting in grading factors rjk employed above, where 

_ r 1 for A; = 1,2,3,4,5,6 , . 

'^'^ ~ 1 n f^v i- — 7 s Q ■ V^j 



^0 for = 7, 8, 9 

In the isospin limit m„ = m^, the quark mass matrix of SU {6\3) is given by 

nig = diag(m„, m„, m^, rrij, rrij, nir, rriu, rriu, rUs). (5) 

Notice that each valence quark mass is degenerate with the corresponding ghost quark's 
mass. This equality maintains an exact cancellation in the path integral between their re- 
spective determinants. Putting only valence quarks in the external states, the contributions 
to observables from disconnected quark loops come from the sea sector of PQQCD. Thus 
one can separate valence and sea contributions as their names suggest, and independently 
vary the masses of the valence and sea sectors. 

The light quark electric charge matrix Q is not uniquely defined in PQQCD [s^. By 
imposing the charge matrix Q to be supertraceless in S'f/(6|3), no new operators involving 
the singlet component are introduced. This can be accomplished with the choice 

Q = diag(g„, g^, g„ qj, qi, g^, qa, g^), (6) 

where to maintain supertracelessness qj + qi + qr = 0. Notice that when the sea quark masses 
are made degenerate with the valence quark masses, QCD is only recovered for the specific 
choice qu = qj = \, and qd = qs = qi = qr -\ — |- As we shall detail below, one can use 
unphysical charges for both valence and sea quarks as a means to determine the LECs. The 
results and problems for the commonly used partially-quenched form of Q proposed in |2^ 
are presented in Appendix O 



B. Mesons of PQxPT 



For massless light quarks, the theory described by the Lagrangian in Eq. (0) has the 
symmetry SU{6\3)l S'[/(6|3)/j (S> f/(l)v, that is assumed to be spontaneously broken to 
SU{6\3)v ® U{l)v in analogy with QCD. One can build an effective low-energy theory of 
PQQCD by perturbing about the physical vacuum state. This theory is PQxPT, and the 
dynamics of the pseudo-Goldstone mesons appearing from chiral symmetry breaking are 
described at leading order in the chiral expansion by the Lagrangian 

£ = :^str (D^S^D^S) + A str (m^S^ + mjs) , (7) 
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where the field 

E = exp(ai^)=^^ (s) 

and the meson fields appear in the SU{6\3) matrix 

The M and M matrices contain bosonic mesons (with quantum numbers of qq pairs and qq 
pairs, respectively), while the x ^i^nd matrices contain fermionic mesons (with quantum 
numbers of qq pairs and qq pairs, respectively). The upper 3x3 block of the matrix 
M contains the familiar pions, kaons, and eta, while the remaining components consist of 
mesons formed from one or two sea quarks, see e.g. ji^. The operation str() in Eq. ((7j) is a 
graded flavor trace. The gauge covariant derivative is deflned by D^T, = d^T, + ieA^[Q, S], 
where is the photon field. To leading order, one finds that mesons with quark content 
qq' are canonically normalized when their masses are given by 



m 



2 



4A 



qq j-2 



— {mq + niqi). (10) 



The fiavor singlet field is $o = str($)/v^, and because PQQCD has a strong axial 
anomaly f/(l)^, the mass of the singlet field has been taken to be on the order of the chiral 
symmetry breaking scale, and thus $o has been integrated out of Eq. ((Tj). The resulting 
rj two-point correlation functions, however, deviate from their familiar form in xPT. In 
calculating the electromagnetic properties of heavy hadrons, the results do not explicitly 
depend on the form of the fiavor- neutral propagator. 

C. Baryons containing heavy quarks in PQ^PT 

Heavy quark effective theory (HQET) and PQxPT can be combined to describe the 
interactions of heavy hadrons and pseudo-Goldstone mesons. Let niq denote the mass of 
the heavy quark, where Q = c oi h. In the limit mq oo, baryons containing a heavy 
quark and two light quarks (and indeed all heavy hadrons) are classified by the spin of their 
light degrees of freedom, s^, because the heavy quark's spin decouples from the system. The 
inclusion of baryons containing a heavy quark into xPT was carried out in 3 2^ ^'^^ 



their magnetic moments were calculated to one-loop order in |4(J,|41|. The quenched chiral 
theory for baryons with a heav y q uark was written down in ||42|] and the partially quenched 
theory was investigated in (jgl-lisj. 

Consider first the baryons with sg = 0. To include these spin-| baryons into PQxPT, we 
use the method of interpolating fields to classify their representations of SU {6\3) [l^ 
23,113143. The Si = baryons are described by the field 



^ As we largely work to leading order in l/mg, the label Q is omitted from all baryon tensors. 
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where i and j are flavor indices; a, b, and c are color indices; and a, /3, and 7 are spin indices. 
The tensor T has the symmetry property 



\niVj'T.. 



(12) 



and forms a 3 9- dimensional representation of SU{Q\'i). These states are conveniently clas- 
sified under the quark sectors acted upon, thus we use the super-algebra terminology 
of [i^li^. Under S'f/(3)vai ® SU{3)sca ® 5'f/(3)ghost, the ground fioor, level A transforms as 
a (3, 1, 1) and contains the familiar = QCD baryon tensor Tij, i.e. T^j 
indices are restricted to the range 1 — 3. With our conventions, we have 



Tij, when the 



1 
71 











(13) 



and the superscript labels the 3-projection of isospin. The first fioor of level A contains 
nine baryons that transform as a (3, 1, 3), while the ground fioor of level B contains nine 
baryons that transform as a (3, 3, 1). The remaining fioors and levels are not necessary for 
our calculation. 

Next we consider the = 1 baryons. These spin | and | baryons are degenerate in the 
heavy quark limit and can be described by one field S!^y The interpolating field for these 
baryons is 



S: 



M:7 



Q 



7,c 



3,b a, a 



£abc{C'j'')aP, 



(14) 



where the tensor satisfies the symmetry 



(15) 



and makes up a 42-dimensional representation of SU{6\3). Under S'f/(3)vai ® SU{3)sea. ® 
SU (3)ghost; the ground fioor, level A transforms as a (6, 1, 1) and contains the familiar S£ = 1 



QCD baryon tensor S^j, i.e. S^j 
The QCD fiavor tensor is 



S^,, when the indices are restricted to the range 1 — 3. 



ij ) 



(16) 



with 



B 



1 vO 
V2^Q 



1 vO 



V2-Q ^ 
I- 



(17) 



and similarly for the B*^. Above is the velocity vector of the heavy hadron, and we 
suppress velocity labels on all heavy hadron fields throughout. Again the superscript on 
these states labels the 3-projection of isospin. The first fioor of level A contains nine baryons 
that transform as a (3, 1, 3), while the ground fioor of level B contains nine baryons that 
transform as a (3, 3, 1). The remaining fioors and levels are not necessary for our calculation. 
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The free Lagrangian for the T and iS^ fields is given by 

£ = -z (S^'v ■ DS^) + A (5^5^) + Ai iS^^M+S^) + A2 (5^5^) strA^+ 

+1 [Tv ■ DT) + A3 (TM+T) + A4 (TT) str7W + . (18) 

We have employed ()-notation for flavor contractions that are invariant under chiral trans- 
formations. The relevant contractions can be found in j2^]. The mass of the T field has been 
absorbed into the static phase of the heavy hadron fields. Thus the leading-order mass split- 
ting A appears as the mass of the S^. This splitting remains finite as {wig — > 00, niq ^0} 
and cannot be removed by field redefinitions due to the interaction of the T and fields. 
We treat A ~ 711^^ in our power counting. The Lagrangian contains the chiral- covariant 
derivative D^, the action of which is identical on T and fields and has the form [29'] 

{D^r\^ = d^%, + vl^%,, + (19) 

The vector and axial-vector meson fields are defined by 

= 1 (^^M^t + ^t^M^) ^ A>^='- [id^i^^ - i'^d^i) . (20) 
The mass operator has the usual definition 

-^+ = ^(em,e + eHO- (21) 

Adding electromagnetic interactions to the covariant derivative can be accomplished by the 
replacement + i^-A^^Qq^ij /'^ + Qij) in Eq- (HH), where Qq is the heavy quark 

charge and the factor of one-half normalizes the action of the Kronecker delta. Gauging of 
the axial- vector fields is not necessary to the order we work. 

The Lagrangian that describes the interactions of the T and fields with the pseudo- 
Goldstone modes is given at leading order by the Lagrangian 

£ = ig2 (S^v^'AfS'') e^,,^ + V2g; [(TA'^^^) + [S^A^T)] . (22) 

The LECs appearing above, the Aj, (72, and (73, all have the same numerical values as those 
used in SU{?>) heavy hadron xPT. This is described in Appendix IbI 

III. ELECTROMAGNETIC PROPERTIES OF THE = BARYONS 

In this section, we calculate the charge radii of the = baryons in PQxPT to leading 
order in the heavy quark expansion. Additionally the magnetic moments are calculated to 
next-to-leading order in the heavy quark expansion in order to ascertain the leading light- 
quark mass dependence. Baryon matrix elements of the electromagnetic current can be 
parametrized in terms of two form factors Fi and F2. In the heavy hadron formalism, such 
a decomposition is 

{T{p')\J^\T{p)) =u{p') 

with q = p' — p a.s the momentum transfer. The Dirac form factor Fi{q^) is normalized at 
zero momentum transfer to the baryon charge Qt in units of the electron charge e, namely 



v^F^{q^) + 



ia^^qy 
2Mt 



F2{q' 



u{p), 



(23) 
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Fi(0) = Qt- The baryon charge is a sum of the heavy quark charge Qq, and the total charge 
of the hght degrees of freedom Qt, i-e. Qt = Qq + Qt- Keep in mind that the hght-quark 
electric charge matrix leaves the valence quarks with arbitrary charges, hence Qt and Qt 
need not have their physical values. Ignoring the Dirac contribution, the magnetic moment 
is given by 

/i = F2(0), (24) 

and the Dirac charge radius by 



,2^ 



(25) 



These electromagnetic observables receive two types of contributions: short range contri- 
butions, which arise from local operators in the effective theory; and long range contribu- 
tions, which arise from loop graphs involving the pseudo-Goldstone modes. The local oper- 
ators characterizing the short distance electromagnetic interactions of heavy quark baryons 
come in two forms. First there are electromagnetic operators corresponding to the current 
interacting with the heavy quark. These operators can be found by matching the HQET 



operators onto the chiral theory |46l . I47l | and thus have fixed coefficients. Secondly there 
are operators which characterize the short distance interaction of the current with the light 
degrees of freedom. Chiral symmetry dictates only the form of these operators, leaving them 
with undetermined LECs. 

The HQET part of the short distance dynamics comes from the anomalous magnetic 
moment and charge radius of the heavy quark. In PQyPT the magnetic moment operator 
is matched in the S£ = baryon sector onto the term 



461 4 



Arar 



{Ta,,T) (26) 



where the heavy quark magnetic moment yUg consists of the Dirac piece plus the anomalous 
part, i.e. = Qq + ^""g^'^'' + . . .. The charge radius contribution from the photon coupling 
to the heavy quark arises from the electromagnetic Darwin term 0, [i^ that appears in the 
HQET Lagrangian as 

^ = -^Qv,D,F^''Q. (27) 
In the Se = baryon sector, this HQET operator matches onto the PQxPT term 



(TT) (28) 



which contributes to the baryon charge radii. 

Now we assess short-distance contributions to electromagnetic observables from the light 
degrees of freedom. In PQxPT the magnetic moment contribution from the light quarks is 
suppressed by AQcn/fng because the light degrees of freedom have se = 0. Thus the leading 
magnetic moment operator has the form 

'^^^^^^ (TQa^^r) F^^ (29) 



AA^niQ 



^ We use F^'' = d^A"^ — d'^A^ for the electromagnetic field-strength tensor. 
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FIG. 1: Loop diagrams contributing to the charge radii of the si = baryons. Mesons are denoted 
by a dashed hne, flavor neutrals (hairpins) by a crossed dashed hne, and the photon by a wiggly 
hne. A thin (thick) sohd hne denotes an = {s£ = 1) baryon. The diagrams in the flrst 
row contribute to the Dirac form factor. The remaining diagrams with a photon have no q^- 
dependence. These, along with the wavefunction renormalization diagrams in the last two rows, 
ensure non-renormalization of the electric charge. 

The local contributions in PQxPT to the charge radii of the = baryons arise from the 
dimension-six operator that is contained in the Lagrangian 

^ = -^{rQr)v^d,F'^''. (30) 

X 

Although the quark charges are arbitrary, the LECs fiT and ct of these two SU {6\3) PQxPT 
operators have the same numerical values as the LECs in three flavor xPT. We make this 
connection in Appendix iBl 

The long distance contributions to these electromagnetic observables arise from the one- 
loop diagrams depicted in Fig. [H The diagrams involve the pseudo-Goldstone mesons and 
the vertices are generated from the leading-order interaction Lagrangian Eq. (j221)- These 
loops, however, do not contribute to the magnetic moments of the S£ = baryons. The 
leading light-quark mass dependence of the magnetic moments enters from heavy quark 
symmetry breaking terms. We shall address these terms below. By contrast, the charge 
radii receive contributions from the one-loop diagrams depicted in the figure. These dia- 
grams are divergent and the resultant scale dependence is absorbed by a counterterm of 
the form in Eq. (jHUj) . In this calculation and throughout this work, we use dimensional 
regularization with a modified minimal subtraction scheme, where we have only subtracted 



terms proportional to 

- -7i5 + l + log47r. (31) 

e 

To compactly write out the g^-dependence of the form factors, we define the abbreviation 




P^=yi-x(l-x)^. (32) 
Combining the local and loop contributions, for the Dirac form factor we find 







X, 

2 



TTl 



2 

J{m^Pcp, A, /i) - J{m^, A, /i) - -g^ 



2 

-x(l - x)g^ g{m^P^, A, /i) 



(33) 



The non-analytic functions arising from loop integrals have been abbreviated in the above 
expression by 



Q(m,6, a) = log— ;r log (34) 

and 

:r(m, 5, /i) = (m^ - 26') log ^ + 26VW^' log " ' "^[^ \ (35) 

The sums in Eq. flS^ are over loop mesons of mass m^. The coefficients appear in 
Table HI and depend on the particular = baryon state T. Notice the charges qj and qi 
do not explicitly appear in these coefficients. In the isospin limit, they always enter in the 
combination q^ + qi which is the same as — by supertracelessness of Q. 

From the table, we see that with efficacious choices for the charges of valence and sea 
quarks contributions from certain loop mesons can be eliminated thereby simplifying the 
chiral extrapolation. For the Aq baryon, we can eliminate all but one loop meson; while 
for the Sq baryons, we can reduce the competition from four loop mesons down to two. 
Alternately there is a choice of sea charges which decreases computation time. One can 
preserve the supertracelessness of the charge matrix with the choice qj = qi = qr = 0. 
For this choice, the lattice practitioner does not calculate closed quark loops with current 
insertion. Furthermore either of the valence charges can be taken to zero, and the LEGS can 
still be extracted from the resulting electromagnetic properties. Once the LECs are known, 
physical predictions can be made. 

For the magnetic moments of the se = baryons, the leading light-quark mass dependence 
arises at next-to-leading order in the heavy quark expansion. Terms contributing to the 
magnetic moments are due to the mass-splitting of the B and B* baryons ^] , as well as 
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TABLE I: The coefficients in SU{6\3) PQxPT. Coefficients are listed for the se = baryon 
states T, and are grouped into contributions from loop mesons with mass m^. 





ju 


ru 


js 


rs 




qu + qd + qr 
qu ~i~ 2'?'' 
qd + \qr 


\{qu + qd) - qr 
\{qu - qr) 
\{qd - qr) 




qs + ^qr 
qs + \qr 




\{qs - qr) 
\{qs - qr) 




FIG. 2: Loop diagrams contributing to the leading light-quark mass dependence of the 5^ = 
baryons' magnetic moments. Mesons are denoted by a dashed line, the photon by a wiggly line, 
and a thin (thick) solid line denotes an S£ = (s^ = 1) baryon. Heavy quark symmetry breaking 
vertices are depicted by squares and arise from terms in Eq. 



from vertices generated from spin-symmetry breaking operators [41]. Both of these effects 
break the U{2Nh) heavy-quark symmetry and are described in PQxPT by the operators 

A2 _ A _ _ 

£ = x^£R [S^^ia^^S") - V2\g^^ [{TA^ia^.S'') + iS^'ia^^A'^r)] , (36) 
mq niQ 

where we have only kept terms relevant for our calculation. The dimensionless parameters 
A and \g have the same numerical values as those in 5't/(3) xPT, see Appendix |B] The 
terms in this Lagrangian lead to new vertices that contribute to the magnetic moments of 
the S£ = baryons. Diagrams involving these heavy quark symmetry breaking vertices are 
shown in Fig. |21 

Aside from these loops, we have additionally contributions from the local operators in 
Eqs. (I26p and (j29p . Thus at next-to-leading order in the heavy quark expansion, we find the 
Pauli form factors of the S£ = baryons to be 

[' dx\g{m^P^,A,fi) + 2\ogP^+l\. (37) 
4vrV ^ Jo ^ J 

The new non-analytic function J^{m, 5, /i) entering above is given by 



J^(m, 5, /i) = -5 log ^ + - m2 log \" "\ ■ i^^) 

/i + — + 

The coefficients are the same as for -Fi(g^) appearing in Table H] and hence the charges 
can be adjusted in the same manner to simplify the chiral extrapolation or to decrease the 
lattice simulation time. 



2 A _ ,/A2 



6 — \/6'^ — m2 + ie 
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IV. ELECTROMAGNETIC PROPERTIES OF THE se = I BARYONS 



In this section, we calculate the charge radii, magnetic moments and electric quadrupole 
moments of the = 1 baryons in PQxPT. Similar to the above decomposition for the se = 
baryons, the current matrix elements of the spin-i B baryons have a decomposition 



{B{p')\J^\B{p))=u{p' 



2M, 



B 



u{p), 



(39) 



where q = p' — p a.s the momentum transfer. The spin-| matrix elements of the electromag- 
netic current J'' for the B* baryons can be parametrized as 



{B:{p')\j'^\B:ip)) = -u,{ji)o^^^uM. 



(40) 



where M/i(p) is a Rarita-Schwinger spin vector for an on-shell heavy baryon. This spin vector 
satisfies the constraints v^Upip) = 'j^u^lp) = 0. The tensor O'^'"^ can be written in terms of 
four independent form factors 
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fJ,U 



vPF*{q^ 



ia'^^qr 



2M, 



B 



F;{q' 



+ 



(2Mb) 



v'Gliq') + 



ia^'^qr 



2M 



G*J 



B 



(41) 



To avoid notational confusion within this section, we have appended a superscript * to 
denote the B* baryon form factors. Extraction of these form factors for the B* baryons 
requires a non-trivial identity for on-shell Rarita-Schwinger spinors [H^l. For our purposes, 
the identity takes the form 



q^ 



2M, 



B 



M 



B 



g"gVU/3(p). (42) 



The conversion of the above covariant vertex form factors to multipole form factors is detailed 
in 113. 

The Dirac form factors Fi{q^) and F^{q^) for the B and B* baryons are normalized at 
zero momentum transfer to the baryon charge in units of the electron charge e, namely 
Fi(0) = F^ifi) = Qb- The baryon charge is a sum of the heavy quark charge Qq, and 
the total charge of the light degrees of freedom Q^, i.e. Qb = Qq + Qs- Keep in mind 
that the light-quark electric charge matrix leaves the valence quarks with arbitrary charges, 
hence Qs and Qs need not have their physical values. The Dirac charge radii for have the 
form < >= 6F[{0), for the B baryons, and < r*^ >= 6-Fj*'(0), for the B* baryons. The 
magnetic moments of the B and B* baryons are defined through their Pauli form factors: 
fi = ^2(0) and /i* = F2{0), respectively. The electric quadrupole moments Q of the B* 
baryons are given by 

Q = -^GUO). (43) 

To the order we work in the heavy quark expansion, the form factor G'^iq'^) = because 
the light degrees of freedom have si = 1. Hence the magnetic octupole moments vanish; see 
Appendix O for further discussion. The relations between the electromagnetic form factors 
of the B and B* baryons are also discussed in Appendix O 

The electromagnetic observables of the B and B* baryons receive short and long range 
contributions in the chiral effective theory. The short range contributions stem from both 
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the electromagnetic interaction of the heavy quark and the electromagnetic interaction of 
the light degrees of freedom. As above, the HQET part arises from the magnetic moment 
of the heavy quark as well as the electromagnetic Darwin term. These contributions are 
encoded in the effective theory through the heavy quark magnetic moment operator 



4m, 



Q 



(5V^,5,) (44) 



where /xg = Qq + ^"^^^q) + . . ., and the heavy quark charge radius operator 

^ = (^"'5.) ^AF"'- (45) 

The electric quadrupole moments of the B* baryons cannot receive any contribution from 
the electromagnetic interactions of the heavy quark. 

In PQxPT the leading contribution from the light degrees of freedom to the magnetic 
moments of the B and B* baryons is contained in the Lagrangian 



Ax 



(S^QS,) F^r (46) 



To be complete to the order we work, the LEG fis must be treated as a linear function of A, 
i.e. /is = fj^g^ + -^f^^s^ ■ Moreover we absorb the finite contributions from loop integrals into 

the parameter Determination of this LEG requires the ability to vary A and for this 
reason we treat the A dependence only implicitly. The magnetic moment coefficient yUs is 
the only LEG that must be treated in this fashion to the order we are working. The charge 
radii receive local contributions from the dimension-six operator contained in the term 

£ = ^(5"Q5,)^;^9,F'^^ (47) 
and finally the electric quadrupole moments arise from 

£ = {S^^QS^^) v^d.F^'^. (48) 

X 



The brackets {. • .} produce the symmetric traceless part of Lorentz tensors, i.e. O^^^^ = 

magnetic octupole moment operator enters at 0{l/mQA^) 
which is beyond the order we work. Although the quark electric charges are arbitrary, the 
electromagnetic LEGs for the = 1 baryons, /i^, cs and Q5, have the same numerical values 
as in SU (3) xPT, see Appendix |B] 

The long distance contributions to the electromagnetic observables of the B and B* 
baryons arise from the one-loop diagrams depicted in Fig. IIVI These diagrams involve 
the pseudo-Goldstone mesons which enter from the vertices generated by the leading-order 
interaction Lagrangian Eq. (j22j) . Galculation of these diagrams along with the tree- level 
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FIG. 3: Loop diagrams contributing to the electromagnetic moments and charge radii of the = 1 
baryons. Mesons are denoted by a dashed line, flavor neutrals (hairpins) by a crossed dashed line, 
and the photon by a wiggly line. A thick (thin) solid line denotes an = 1 (s^ = 0) baryon. 
The diagrams in the first row contribute to the electromagnetic moments and charge radii. The 
remaining diagrams with a photon have no g^-dependence. These, along with the wavefunction 
renormalization diagrams in the last two rows, ensure non-renormalization of the electric charge. 

TABLE IL The coefHcients in SU{<o\^) PQxPT. CoefHcients are listed for the = 1 baryon 
states B, and are grouped into contributions from loop mesons with mass m<^. 



JU ru 



rs 




'^'Qu + qr Qu- Qr 

qu + qd + Qr ^{qu + qd) - qr 
'^qd + qr qd - qr 



\{qu - qr) 

\{qd - qr) 




gs + 2 9r 2^qs-qr) 

qs + \qr \{qs-qr) 
+ qr qs- qr 



























14 



contributions for the B baryons yields 



Qb+ Qq- 



-y 

2 f2 



6 



(4vr/)' 



log Pa, 



2 



9l 



dx 



J{m^P^, -A, /i) - J{m^, -A, /x) + -g^ 

y 



- - G{m^P^, -A, ^u) 



for the Dirac form factor and 

Mb 



2^ Mb AMb^ , ^IMb ^ b f\ p 

F2(g) = — /^Q + /^^^2- + ^2.«. X ' ' 



B 



A, Jo 



(49) 



(50) 



for the Pauh form factor. The function J-'{m,S,n) appearing above is given in Eq. (jHHj) . In 
the above expressions, the sums run over loop mesons (p with mass m^. The coefficients 
are listed in Table [H] and depend on the particular baryon state B. Notice the charges Qj 
and qi do not explicitly appear in these coefficients. In the isospin limit, they always enter 
in the combination qj + qi which is the same as — g,. by supertracelessness of Q. 
Carrying out the calculation for the B* baryons, we obtain 



F*{q'' 



SttV 



-y 

2 f2 



dx 



q 



6 



dx|2[m5-2x(l-x)g2]logP^-ig2 



mh 
log^ + 1 



+ 



9l 



E 



(4vr/)2 ^ 
for the Dirac form factor 



a 



B 



dx 



J{m^P^, -A, /i) - J{m^, -A, /i) 



p*. 2x Mb ^ 2Mb^ qIMb^ b \ p 
"^Q 167r/^ ^ Jo 



glM. 



B 



y^^ dxJ='{m^P^,-A,jj,) 
J. Jo 



(51) 



(52) 
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for the Pauli form factor, and finally 



Mb 



A. 



M 



B 



A 



X 



B 



1 

mi 



ai / dx 





x(l - x) g{m^P^, -A, yu) + - 

oJ 



(53) 



The coefficients appearing in the expressions for the B* baryon form factors are listed 
in Table ITTl and are identical for baryon states B and B*. Taking the rriQ — > oo limit, the 
relations between the B and B* electromagnetic form factors (derived in Appendix [CJ are 
satisfied by our one-loop PQxPT results. 

As with the se = baryons. Table ITTl shows that the chiral extrapolations can be simplified 
with efficacious choices for the charges of valence and sea quarks. For the Eq and Qq 
baryons, we can eliminate all but one loop meson; while for the baryons, we can reduce 
the contributions from four loop mesons down to two. In SU {6\3) one can again decrease 
the computation time by choosing the sea quarks to have vanishing electric charges as this 
preserves the supertracelessness of the charge matrix. Thus all closed quark loops with a 
current insertion vanish. Furthermore for the Eg, and Sg baryons, either of the contributing 
valence charges can be taken to zero, and the LEGS can still be extracted from the lattice 
data. For the Eg^, Eg^, and fig baryons, this choice is not possible since these states 
have two valence quarks of the same fiavor. Nonetheless, LECs extracted from calculations 
involving one baryon can be used to test PQxPT predictions for the remaining baryons. 



V. SUMMARY 



Above we have calculated the electromagnetic properties of baryons containing a heavy 
quark. For baryons with light degrees of freedom se = 0, we obtained the charge radii at NLO 
in the chiral expansion. The magnetic moments of these baryons vanish at leading order 
in the heavy-quark expansion and thus to obtain the leading light-quark mass dependence 
we worked to NLO in the heavy quark expansion. For the = 1 baryons, we obtained the 
charge, magnetic moment and electric quadrupole moment form factors at NLO in the chiral 
expansion. Expressions for these quantities are derived in PQxPT for three light fiavors in 
the main text, and two light fiavors in Appendix^ Additionally the corresponding xPT 
results for three and two light fiavors are presented in Appendix El Knowledge of the quark 
mass dependence of these observables is essential to extrapolate lattice QCD data to the 
physical light quark masses. Additionally the effective field theory predicts the momentum 
transfer dependence of the electromagnetic form factors near zero recoil. This information 
should be utilized to perform the zero-momentum extrapolation of lattice data in order 
to extract static electromagnetic properties. Furthermore there are relations between the 
various sextet baryon electromagnetic form factors that follow from the decoupling of the 
heavy and light spins; these are found in Appendix O 

There are a number of interesting points that arise when considering the electromagnetic 
properties of baryons containing a heavy quark as compared to baryons containing only 
light quarks. Firstly the choice of the heavy quark lattice action allows for more freedom in 
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separating contributions to baryon properties. For the octet and decuplet baryons formed 
from three hght quarks, in order to have a consistent power counting the average inverse 
baryon mass 1/M becomes an expansion parameter 51, S^l- On the lattice, however, there 



is no direct control over this expansion, whereas the choice of the HQET action leads to 
control over the l/mg expansion. This in turn can provide information about the chiral 
regime. Consider the S£ = baryons. The form factors Fi{q^) and F2{q^) enter at different 
orders in the heavy-quark expansion. Thus if the light-quark mass dependence from xPT 
can be trusted, then increasing the heavy quark mass should diminish the variation of -^2(5^) 
and isolate the leading loop contributions to Fi(g^). If this behavior is not observed, one is 
certainly not in the chiral regime. 

A further contrast to the baryons with only light quarks, is the case of the Aq and TPq 
baryons. In SU (2) flavor, the loop contributions to the electromagnetic properties of these 
baryons vanish due to cancellations between 7r+ and vr" loops. Such complete cancellations 
do not occur for the nucleons and deltas.^ From the point of view of chiral physics, this 
is perhaps uninteresting as there is no quark mass dependence. This situation is true at 
NLO in the heavy quark expansion as well, because the flavor structure of the heavy quark 
symmetry breaking operators is identical to the LO operators. If there is observed quark 
mass variation of the electromagnetic properties of the Aq and Eg, then one is not in the 
chiral regime. Even though there is no quark-mass dependence in QCD, the couplings g2 and 

can still be extracted from considering the Aq and TPq in PQ^PT, because the arbitrary 
quark electric charges can be chosen so that loop diagrams do contribute. 

A feature of three-flavor electromagnetism, is the tracelessness of the charge matrix. This 
carries over into the partially quenched theory as supertracelessness, and enables the LECs 
to be extracted from lattice data in which the sea quarks have vanishing charges. This 
reduces the computation time, as closed quark loops with current insertion vanish. We have 
accordingly used a charge matrix which clearly separates valence and sea contributions; see 
Appendix ^ Similar electroweak operators should be utilized for other observables. 

Lattice QCD calculations will enable an understanding of the electromagnetic properties 
of hadrons in terms of quark and gluon degrees of freedom. These observables paint an 
intuitive picture of the low-energy electromagnetic structure of hadrons. Foreseeable lattice 
calculations of these properties, however, will require extrapolation in the light quark masses. 
The expressions derived in this paper for heavy hadrons in xPT and PQxPT allow such 
extrapolations, as well as the zero-momentum extrapolation. The PQxPT expressions are 
not complicated by new non-analytic functions arising from the flavor-neutral propagator. 
Moreover a study of the electromagnetic properties of baryons containing a heavy quark 
allows the exploration of our ability to extrapolate lattice QCD on two fronts: the light-quark 
regime and the heavy-quark regime. Control over the latter scale is not a feature for nucleons 
and deltas treated as heavy baryons. Lastly the exploitation of arbitrary electroweak quark 
charges allows for calculational simplifications to extract physical low-energy constants. 



^ There is the case of the spin-| T,*-^ baryon in SU{i) for which the pion and kaon loops individuaUy sum 
to zero. The local electromagnetic operators, however, are all proportional to the baryon charge, thus to 
one-loop order all electromagnetic observables of the vanish. 
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APPENDIX A: ELECTROMAGNETIC PROPERTIES IN SU{A\2) 

Heavy hadron chiral perturbation theory in the baryon sector for SU{A\2) parallels that 
of SU {6\3) in Sec. IHI The partially quenched heavy hadron Lagrangian has been written 
down in [29]. In this appendix, we briefly review the setup of S'f/(4|2) PQQCD and PQxPT. 
We then present the calculation of heavy hadron electromagnetic properties for this graded 
flavor group. 



1. PQxPT for 5^7(412) 

In SU{4:\2) PQQCD, the light-quark sector Lagrangian appears as 

6 

£= 5^ g^.(z^-m,).,gfe, (Al) 
j,k=i 

and differs from the usual QCD Lagrangian by the inclusion of four extra quarks; two bosonic 
ghost quarks, (m, d) , and two fermionic sea quarks, (j, /) , in addition to the light physical 
quarks {u, d). The six quark fields transform in the fundamental representation of S't/(4|2). 
They appear in the six-component vector 

q = {u,d,j,l,u,d)^. (A2) 

The quark fields obey the graded equal-time commutation relations in Eq. (jH)) but with the 
grading factor rjk now defined by 

' \ for A; = 5, 6 ^ ' 

The SU{4\2) mass matrix with non-degenerate quarks is given by 

ruq = diag(m„, m^, mj, mi, ruu, ma), (A4) 

where the ghost quarks remain degenerate with their valence partners. Defining ghost and 
sea quark charges is constrained only be the restriction that QCD be recovered in the limit 
of appropriately degenerate quark masses. The general form of the charge matrix we choose 
is 

Q = diag {qu,qd,qj,qi,qu,qd), (A5) 

which is not supertraceless. In the limit rrij —>■ rriu and mi — > m^, QCD with two light flavors 
is recovered only for the values qu = <ij = \, and qd = qi = —\- The results for and problems 
with the commonly used partially-quenched charge matrix are presented in Appendix iDl 
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For massless light quarks, the Lagrangian in Eq. ()A1|) has a graded symmetry SU (4|2) l ® 
SU{4:\2)ji^U{l)v, which is assumed to be spontaneously broken to SU{4:\2)v ®U{l)v- The 
Lagrangian describing the pseudo-Goldstone mesons of this theory is identical in form to 
Eq. ((7j). The meson fields appearing in Eq. (jU)) are replaced by an SU{4\2) matrix $. The 
upper 2x2 block of the matrix M contains the familiar pions, see e.g. 0|. The flavor 
singlet field that appears in ^[/(4|2) is defined to be % = str($)/v/2. As before, the mass 
of the singlet field mo can be taken to be on the order of the chiral symmetry breaking scale. 
In this limit, the 7] two-point correlation functions deviate from their familiar form in xPT, 
and their explicit forms are not needed for expressing our results. 

To include baryons with one heavy quark into the theory, we use the same interpolating 
fields, Eqs. and (fT^ . The S£ = baryons are then described by the field T^j which forms 
a 1 7-dimensional representation of S't/(4|2). The baryon tensor of QCD Tij is contained as 
Tij = Tij, when the indices are restricted to the range 1 — 2. In our conventions, we have 



T = — ( ^ Ac 
V2\-Aq 



(A6) 



The Sc = I baryons are described by S^j which makes up a 19-dimensional representation 



of SU{A\2). The baryon tensor of QCD S^^j is embedded as S^^j = S^^, when the indices are 
restricted to the range 1 — 2. Here the QCD flavor tensor S^^ is given in terms of Bij and 
B*j as in Eq. dTHI), but with 




and similarly for 3*^". The superscript on these states labels the S-projection of isospin. The 

remaining states relevant to our calculation have been classified in j29l]. 

The free Lagrangian for the T and fields is the same as Eq. fll8|) and the Lagrangian 

that describes the interactions of these fields with the pseudo-Goldstone modes is given 

by Eq. (jUJ. The LECs appearing in the Lagrangian, Ai, A2, g2, and (73, all have the 

same numerical values as those used in SU (2) heavy hadron xPT, see Appendix |B] These 

are of course numerically different than those in S'f/(3). The parameters Ag^*^"* and A4^'^^ 

are different and can be related to the unquenched coefficient A4 by matching, namely 
\ _ 1 \LPQ) , \{PQ) 

A4 — TjAq A, 



2. Baryon electromagnetic properties 

Calculation of the baryon electromagnetic properties in SU{4\2) is similar to that in 
SU{Q\?)) because the Lagrangian has the same structure. There is, however, one crucial 
difference because the charge matrix is not supertraceless and hence there are new operators 
contributing at tree level to each of the electromagnetic properties. 

For the Dirac form factor of the = baryons, there is an additional operator that 
contributes at tree level, namely 

£ = (^^) vAF^"" strQ. (A8) 

X 
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TABLE III: The coefficients for the Aq in SU{A\2) PQ^PT. Coefficients are grouped into 
contributions from loop mesons with mass m^. 





ju lu 


jd Id 






Uid-qj) Uid-qi) 



The contribution of this operator to the Dirac form factor of the Aq baryon is 

m{q') = (A9) 

X 

where qji = qj + qi is the supertrace of the charge matrix. The remainder of the form factor 
is given in Eq. (jH^ . however the coefficients are now hsted for the Aq baryon for 5'f/(4|2) 
in Table IIIII For the magnetic moment of the Aq , we have the additional operator 

£ = -^Mq^ (Ta^^r) F^^ strQ, (AlO) 

which leads to a contribution to the Pauli form factor 

5F,(g) = g,,/.^-^, (All) 

and the remainder of the Pauli form factor is given in Eq. (jHTf) with the coefficients for 
SU{4\2). 

From Table IIIH we see that the charges of the valence and sea quarks can be chosen so 
that there are no loop contributions at this order. This dramatically simplifies the partially 
quenched chiral extrapolation and can be anticipated from the SU{2) result, see Appendix iBl 
One can choose Q to be supertraceless to sort out ct from ct, and /xr from Jirp. The specific 
choice qj = qi = is possible, but the result is not sensitive to all the LECs. Thus in SU{4:\2) 
quark disconnected contributions with current insertion must be dealt with to ascertain all 
LECs and hence make physical predictions. 

Considering next the = 1 baryon electromagnetic properties, we must add similar types 
of operators involving the supertrace of the charge matrix. These operators are 

£ = !^ (S^S,) F'^-' strQ + ^ v^d^F'^" strQ - ^ (5^^5^>) v^d^F,'^ strQ. 

A^ 

(A12) 

For the spin-| B baryons, the contributions to the form factors from the above terms are 



XT? f 2\ f . Qs\ 2 



X X' 

5F,{q') = ^q.ijls- (A13) 

The remaining parts of the Dirac and Pauli form factors are the same as in Eqs. and 
flM?|) . respectively. The coefficients for SU{A\2) are listed in Table HVl for the Eg baryons. 
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TABLE IV: The coefficients in 5C/(4|2) PQxPT. Coefficients are listed for the Si 
states B, and are grouped into contributions from loop mesons with mass m^. 



1 baryon 





ju 


lu 


3d 


Id 


y + l 
y-1 


Qu - Qj 

liln - qj) 



Qu - qi 
hiQu - qi) 





\{qd-qj) 

qd - qj 




\{qd - qi) 
qd - qi 



Lastly, for the spin-| B* baryons, the contributions to the electromagnetic form factors from 
the above operators are 



Cs_ 

X 



2A^Mb 



2AI 



2M, 



B _ 



Ax 
4Mb 



-Qji t^s 



/^5 + 



Mb 



Ax V"' A, 



(A14) 



with the remaining parts of these three form factors given in Eqs. ()51|) . and (j^Hj) with 
the coefficients listed in Table HV1 

From Table ITvl we see that for each Eg baryon the charges of the valence and sea quarks 
can be chosen so that there are no loop contributions at this order. Consequently there 
is a simple partially-quenched chiral extrapolation. Compared with the SU{2) result, see 
Appendix El only the Eg has this simplification. One can choose Q to be supertraceless 
to sort out cs from C5, etc. The specific choice qj = = is possible, but again does not 
lead to results that depend on all LECs. Quark disconnected contributions with current 
insertion must be calculated in SU{4\2) to ascertain the relevant LECs and make physical 
predictions. 



APPENDIX B: ELECTROMAGNETIC PROPERTIES IN SU{3) AND SU{2) 
For completeness, we include calculations of the baryon electro mag netic properties in 



xPT. The magnetic moments were calculated to one-loop order in |40|,|4H; the remammg 
properties have not been calculated before. To compare with the main text, we consider the 
case of SU{3) flavor in the isospin limit and to compare with the results of Appendix [XJ we 
consider non-degenerate SU(2). For these calculations, we retain the tensors T^- and S^^ in 
xPT. For the case of SU{?)) flavor the tensors are given in Eqs. (fT^ and (fTTj) . respectively; 
while for SU{2) flavor, they are given by Eqs. ()A6|) and ()A7|) . respectively. 

The free Lagrangian for the Tij and S^j flelds in SU (3) and SU{2) both have the form 

£ = -I (Sf'v ■ DS^) + A (5^5^) + Ai (S^'M+S^) + A2 (S^S^) tiM+ 

+t {Tv ■ DT) + A3 (TA^+T) + A4 (TT) trA^+, (Bl) 

with the exception that A3 = in the two-flavor theory to avoid redundancy 4^. Similarly 
the interaction Lagrangian for Tij and S^^ has the same form for both flavor groups, namely 

£ = ig2 (S^v^A^S'^) E^,p, + 72^3 [(TA^'S^) + (S^A.T)] (B2) 
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TABLE V: The coefficients in SU (3) xPT. Coefficients are listed for the S£ = baryon states 
B, and are grouped into contributions from loop mesons with mass m^. 





TT 


K 







1 

2 




1 

2 





" 2 


1 


1 




2 


2 



The factor a/2 is chosen so that g^i agrees with the Uterature, e.g. [53]. The numerical values 
of the Aj, g2, and are of course different for the SU{?>) and SU{2) theories. 



1. 5C/(3) 

In contrast to the partially quenched theories considered above, the light-quark mass 
matrix in the isospin limit of SU (3) flavor is given by 

mq = diag(m„, m„, m^), (B3) 

and the light-quark charge matrix is 

The local operators contributing to the electromagnetic properties of the = baryons are 
contained in the Lagrangian 

^ = - J^ {Ta,,T) F^'^ - 1% {TT) v.d.F'^^ (B5) 
for the heavy quark contributions and 

^ = _ ^^T^QCD (TQa,,T) F^'' - ^ (TQT) v,d,F^^, (B6) 

for the brown muck. The LECs and ct in 5'f/(6|3) have the same numerical values as in 
S'f/(3) xPT. The heavy quark symmetry breaking operators have the form 

£ = i^^ia^^S'') - V2Xg^^ [{TAHa^^S") + (S^ta^^A'T)] (B7) 

for the 5'f/(3) flavor group. We trivially find the LECs above, A and Ag, have the same 
values as in S'f/(6|3). The Dirac and Pauli form factors can now be determined in SU{?>) 
and they have a form identical to Eqs. (jH^ and (j?fj) . The coefficients are listed for 5'f/(3) 
in Table El 

The local operators contributing to the electromagnetic properties of the si = 1 baryons 
are contained in the Lagrangian 

£ = ^ F^'' + 1% (S^S^) (B8) 
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TABLE VI: The coefficients in SU{3) xPT. Coefficients are listed for the S£ = 1 baryon states 
B, and are grouped into contributions from loop mesons with mass m^. 







K 
i\ 


y + l 


1 


1 







1 
2 




-1 







1 

2 







1 


1 




2 


2 







-1 



for the heavy quark contributions and 

for the brown muck. The LECs ixsi cs, and Qs in SU{6\3) have the same numerical values 
as those in SU{3) xPT. The three form factors Fi[q^), and -^2(5'^) can now be determined 
in SU{?>) for the B baryons and they have a form identical to Eqs. and (j^Uj) . The form 
factors F^{q^), F^i^q^) and G\{q^) for the B* baryons have the same form as Eqs. (jSTJ), 
and ()53p . For each of the form factors for the B and B* baryons, the coefficients are 
listed for 5^7(3) in Table ED 

2. SU{2) 

The light-quark mass matrix in SU{2) flavor is given by 

niq = diag(m„,md), (BIO) 



and the light-quark charge matrix is 



Q = diagQ,-i^. (Bll) 

The operators contributing to baryon electromagnetic properties have the same form as in 
SU{3). One must keep in mind, however, that the LECs appearing in Eqs. fjB6|) . ()B7|) . and 
flB9|) have different numerical values in SU{2). These values are the same as in SU{A\2). 
Additionally there are local electromagnetic operators in SU{2) that involve the trace of the 
charge matrix. These operators have the form 

o _ ^^^^ ^^g^^,, _ e-JlTAQCD ^^^^^^^ 



(^S^SJ) F^'^tiQ + ^ {S^S:) v.d^F^'^tiQ - ^ {S^'^S''^) v^d.F^^tiQ. 
^^x ^\ ^\ 



(B12) 
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The LECs appearing in 5'?7(4|2) are all numerically equal to those in the above Lagrangian 
for SU{2). Including contributions from these operators, the electromagnetic form factors 
of baryons in SU{2) have exactly the same form as those in SU{4\2) in Appendix 1X1 with 

Qji = |. For the Ag baryon, however, the are no net loop contributions since at''^ = 0. For 
the Sq baryons, the value of the coefficient a is given by the isospin: a-n-'^ = I. 

APPENDIX C: FORM FACTOR RELATIONS IN THE HEAVY QUARK LIMIT 

In this Appendix, we derive relations between the electromagnetic form factors and tran- 
sition form factors of the B and B* baryons. These relations follow from the decoupling 
of the spins of heavy and light degrees of freedom in the heavy quark limit. We begin for 
simplicity with the case of the = baryons. 

The electromagnetic current can be written in terms of light quark and heavy quark 
contributions, 

i i 

where the first sum is over the light fiavors while the second is over the heavy fiavors. Written 
this way, our arguments generalize to partially quenched theories. Notice that only the heavy 
quark charge operator survives the heavy quark limit, thus the action of between heavy 
quark states is proportional to unity. Consider the current matrix elements of the = 
baryons. As the spins of the heavy and light degrees of freedom are decoupled, the baryon 
states are direct products of = states and heavy quark states. The action of the current 
operator is always unity in the heavy quark subspace, thus we have 

{T{p')\J^\T{p)) = {Q\l\Q)v^f{q'). (C2) 

Comparing the above form with the current matrix element decomposition in Eq. ()23|) . we 
find Fi(g2) oc f{q^) and F2{q^) = 0. 

Now we consider the current matrix elements of the sextet baryons. In the heavy quark 
limit, the B and B* baryons are degenerate and are described by the = 1 field defined 
in Eq. Schematically current matrix elements of the S field can be written as 

(S,{p')\J''\S,{p)) = ((s, = l,p'\ ® {Q\)r{\Q) ® \se = l,p)), (C3) 

because of the decoupling of heavy and light spins. As the action of in the heavy quark 
sector is unity, current matrix elements of the S-field can be parametrized by 

{S,{p')\r\S,{p)) = -e,{p')0^'"'e,{p), (C4) 

where is the polarization spinor of the S-field that is proportional to the direct product: 
u{v) ®e^{p), where u{v) is the heavy quark spinor and e^{p) is the polarization vector of the 
light degrees of freedom. The polarization spinor satisfies the constraints p ■ e{p) = ^{p) = 
V ■ e{p) = 0. The tensor O'^^'^ has the decomposition 

O.P. = g'^-'v^Mq') + ^iq'g'" - gV^)/2(g') + WKfT2^'9iiq'). (C5) 
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Notice that there are only three form factors for the S'-field because the hght degrees of 
freedom have spin one. The above current matrix element encodes the electromagnetic 
form factors of both the B and B* baryons, as well as their electromagnetic transition form 
factors. 

Comparing the current matrix element of the B* fields in Eq. (j41|) with that of Eq. ()C4|) 
allows us to determine 

Fliq') = h{q% 
F;{q') = /2(g^), 
Gliq') = gi{q'), and 

Gliq') = 0, (C6) 

in the heavy quark limit. Notice that the magnetic octupole form factor 6*2(0'^) vanishes by 
heavy quark symmetry. Carrying out the procedure for the current matrix element of the 
B fields in Eq. we find 

Fi{q') = fi{q') + ^9i{q% and 

F2{q') = ^/2(g'), (C7) 

in the heavy quark limit. Thus we have relations between the B and B* electromagnetic 
form factors, namely 

Fi{q')-F*{q') = ^Gl{q'), (C8) 

and 

F2{q')-lF;iq')=0. (C9) 

These in turn yield relations between the static electromagnetic properties of the B and B* 
baryons 

< r- > - < >= ^, (CIO) 

and 

/i=^/i*. (Cll) 

These expressions are satisfied in the heavy quark limit by our one-loop PQxPT and xPT 
results. 

Lastly we consider the electromagnetic transition matrix element 

(B;{p')\r\B{p)) = u,{p')o'"'^Mp), (ci2) 

where the tensor O^'^ is given by 
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TABLE VII: The coefficients a J in 5C/(6|3) PQxPT for the charge matrix Q in Eq. (jDTJ). Co- 
efficients are hsted for the si = baryon states T, and are grouped into contributions from loop 
mesons with mass m^. The abbreviation qji is defined to be qji = qj + qi. 

TT K Tjs ju ru js rs 



A, 



Q 

2 

Q 



'6 + 2^i' 2 ('^'i' + ^'") 6 + 2'?'^ 3 2^jl 3 2^?^ 3 2^i' 

1 J- Irr J- n i J- i/i _i _ l/i _i _ Irr _1 _ l/i 

'6 '2^1 ''ir) 6 ' 3 2'ii' 6 3 2^^' 



.1 _ la 
6 2y'' 

.1 _ la 

6 2'i'' 



Matching this decomposition onto that in Eq. ()C4|) for the S'-field yields 

G2{q^) = -^^i(g^), and 

Gsiq') = 0. (CM) 
In particular these relations imply that the magnetic dipole transition moment is given by 

Gi(0) = -^/i*, (C15) 

while the electric quadrupole transition moment is given by 

G2(0) = -^Q, (C16) 

and the Coulomb quadrupole moment 6*3(0) vanishes. 

APPENDIX D: QUARK CHARGES 

In this Appendix, we provide results for the more commonly used choice of the charge 
matrix Q. The choice used in the main text has advantages over this commonly used form 
as we explain below. Insertion of different quark charges results only in modification of the 
charge dependent factors in our above results. Thus the only factors that are altered for a 
different choice of Q are the baryon charge Qt and Qb, the charge of the light degrees of 
freedom Qt and Qs, and the coefficients and that depend on the charge of the loop 
meson 0. 

The commonly used form of the charge matrix Q in SU{6\3) PQQCD is ji^ 

Q = diag - J' ■ (^^1) 

When the sea and valence quark masses are made degenerate, one rather elegantly recovers 
QCD independent of the values of the charges qj, qi, and g^- Implementation of this choice 
for the charge matrix, however, is not the most economical in terms of computation time. 
This is because closed valence quark loops with a photon insertion are not canceled by 
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TABLE VIII: The coefficients a| in 5C/(6|3) PQxPT for the charge matrix Q in Eq. JdTJ). Co- 
efficients are Hsted for the S£ = 1 baryon states B, and are grouped into contributions from loop 
mesons with mass m^. The abbreviation qji is defined to be qji = qj + qi. 
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TABLE IX: The coefficients for the Aq in 5C/(4|2) PQxPT for the charge matrix Q in Eq. (lD2l) . 
Coefficients are grouped into contributions from loop mesons with mass m^. The abbreviation qji 
is defined to be qj + qi . 
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the corresponding ghost quark loops, as their charges are not identical. Thus the lattice 
practitioner must calculate closed valence quark loops with photon insertion (using the 
effective charges = | — Q'j, qd = ~\~'ih and = — | — to mimic the partial cancellation 
from ghost quark loops) in addition to closed sea quark loops and closed sea quark loops 
with photon insertion. Closed valence quark loops without photon insertion are of course 
canceled by closed ghost loops because their charges never enter. As valence quarks now 
appear in the closed loops with an operator insertion, the benefits of using lighter valence 
masses than sea masses cannot be maximally obtained. In the foreseeable future, the choice 
of Q above is not ideal for performing partially quenched simulations. 

In order that the valence and sea sectors be separated as their names suggest, one must 
have the ghost charges equal to their valence counterparts. With this choice, closed valence 
quark loops with operator insertion are always canceled by the corresponding ghost loop. 
Thus the only quark disconnected contributions arise from the sea sector and these quarks 
efficaciously can be given larges masses to reduce computation time. Thus a more practical 
choice of the SU{Q\?)) charge matrix is that employed in the main text 

Q = diag(g„, g^, g„ qj, qu gr, gn, gd, g^), 



TABLE X: The coefficients in 5C/(4|2) PQ^PT for the charge matrix Q in Eq. Coefficients 
are hsted for the = 1 baryon states B, and are grouped into contributions from loop mesons 
with mass m^. The abbreviation qji is defined to be qji = qj + qi. 
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where to maintain supertracelessness qj + qi + qr = 0. Notice that when the sea quark masses 
are made degenerate with the valence quark masses, QCD is only recovered for the specific 
choice qu = Qj = I, and qd = qs = qi = qr -\ — |- One can use unphysical charges for both 
valence and sea quarks to determine the electromagnetic LECs. 

To contrast with the calculation in SU{Q\?)) in the main text, we calculate the baryon 
electromagnetic properties using the charge matrix in Eq. ()D1|) . The results have the same 
functional form as in the main text. One must be careful, however, to use the appropriate 
values of the baryon charge, the charge of the light degrees of freedom, and replace the 
5'C/(6|3) coefficients and with those listed in Tables IVTTl and IVIIII Incomplete can- 
cellations between valence and ghost sectors imply that valence- valence meson loops are not 
completely canceled by valence-ghost meson loops. Not surprisingly then these results have 
more complicated chiral extrapolation formula because more loop mesons contribute and 
one can try to simplify the formula by using clever choices for the unfixed charges. 

As with the case of SU{Q\?>), the choice of charge matrix in for S'f/(4|2) 

Q = diag Q,-i,gj,gz,gj,gz^ (D2) 

is computationally intensive. A choice that more readily can be put to use on current lattices 
is 

Q = diag(g„, g^, qj, qi, q^, qd), 

which we used in Appendix ^ This choice mandates that the only closed quark loops 
are those with sea quarks. Carrying out the calculation using Q in Eq. ()D2|) . we find the 
modified 5'f/(4|2) coefficients and listed in Tables HXl and IXf respectively. Again there 
are more contributing loop mesons due to incomplete cancellations between the valence and 
ghost sectors of the theory. 
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